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LECTURE 14. STABILITY 



The notion of stability. Roughly speaking, a system is called stable if its long-term behavior does 
not depend on significantly the initial conditions. 

An important result of mechanics is that a system of masses attached in (damped or undamped) 
springs is stable. A similar result is in network theory. In these notes, we study the differential 
equation of the form 

(14.1) y"+py' + gy = f^t), 

where p, q are constatnts and f{t) represents the external forces. 
We learned that the general solution of (14.1) has the form 

(14.2) y = civi + C22/2 + Vp, 

where ci, C2 are arbitrary constants and yp is a particular solution of (14.1); ciyi + C2y2 is the 
complementary solution, that is, the general solutions of the homogeneous equation (14.1) with 
fit) = 0. 

The initial conditions determine the values of ci and C2. Thus, we say the system (14.1) is stable 
if ciyi + 022/2 ^ as t ^ cx) for any coice of ci and C2. 

If (14.1) is stable then yp is called the steady-state solution and cij/i + C2?/2 is called transient. 
Physically, in a stable system, the output is the sum of a transient term, which depends on the 
initial conditions, but whose effects die out over time, and a steady-state, which represents the 
response of the system to the input f{t) after a long time. 

Stability conditions. We study under what circumstances the differential equation Ly = /, where 

(14.3) L = + piD""^ + +p^_^D+pn, 

where pj are constants, is stable. 

Definition 14.1. The differential equation Ly = f, where L is given in (14.3) is called: 

(i) asymptotically stable if every solution of Ly = tends to zero as t ^ oo; 

(ii) stable if every solution of Ly = remains bounded a t — > oo; 

(iii) unstable if it is not stable. 

We note that stability concerns only the behavior of the solutions of the corresponding homo- 
geneous equation Ly = 0. 

When f{t) = 0, then a steady-state solution is y = 0. In this case, the system is stable if small 
initial departures from the steady-state remain small with the lapse of time. 

By definition, Ly = f is asymptotically stable if every basis solution of Ly = tends to zero 
as t ^ oo and it is stable if the basis solutions remain bounded. In view of the characteristic 
polynomial of L and the fundamental theorem of algebra, we write 

L = {D- Ai)'=H^ - ^2)^' ■■■{D- 
where \j G C are all distinct and ki + k2 + + km = n. 

Exercise. The general solution of the homogeneous equation Ly = is given by 

yit) = ci(t)e^^* + C2(i)e^2* + c„(t)e^'"*, 
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where Cj{t) is an arbritary polynomial of degree kj — I. 

Exercise. If r is a nonnegative integer and A G C, show that 

lim 1^6^*1 =0 if ReA < 0. 

Therefore, Ly = f is asymptotically stable if ReAj < for all j, and it is stable if ReXj < or 
ReAj = and kj = 1. 

We summarize the result. 

Theorem 14.2. The differential equation Ly = f is asymptotically stable if every root of the characteristic 
polynomial of L has a negative real part, and it is stable if every multiple root has a negative real part and 
no simple root has a positive real part. 

Example 14.3. We consider the second-order differential equation 

(14.4) y" + py' + qy = 0, p,q are constants. 

We recall that the discriminant A = — 4g tells us about the nature of the solutions, and hence 
about the stability of (14.4) 

Ifq < then A > and the characteristic polynomial A^ +pX+qhas two real roots with opposite 
signs. Therefore, (14.4) is unstable. 

If p < then at least one root of the characteristic polynomial must have a positive real part. 
Hence, (14.4) is unstable. 

If p = and q > 0, then (14.4) reduces to y" + qy = with q > 0. Hence, it is stable but 
asymptotically stable. 

Finally, let p > and q > 0. If A ^ then the roots of the characteristic polynomial have 
negative real parts, and (14.4) is asymptotically stable. If A > then A = p"^ — 4q < p"^ and thus 
a/a < p. Therefore, (14.4) is asymptotically stable. 

In summary, (14.4) is asymptotically stable if and only if p > and q > 0, and stable if and only 
if p ^ and q > 0. 

Stability of higher-order differential equations. The above example phrases the stability crite- 
rion for (14.4) in terms of the coefficients of the equation. This is convenient since it does not 
require one to calculate the roots of the characteristic polynomial. 
For higher-order equations, 

(14.5) + piy^""^^ + • • • + Pn-iv' + PnV = 0, Pj are constants, 

it is not too hard to show that if (14.5) is asymptotically stable then pj > for all j (Exercise). 
But, the converse is not true (Exercise). For the implication of a criterion for coefficients of (14.5) 
for stability, the coefficients must satisfy a more complicated set of inequalities, which we state 
without proof in the following. 

Routh-Hurwitz Criterion for Stability. The differential equation (14.5) is asymptotically stable if 
and only if in the determinant 

pi 1 ... 

P3 P2 Pi 1 ••• 

P2n-1 P2n~2 Pn 
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where pk = if k > n, all of its n principal minors, that is, the subdeterminar\ts in the upper left 
comer having sizes respectively 1,2, ... ,n, 



Pi, 
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are positive. 



Exercise. We consider 



{D'^ + 2D^ + 6L>2 + 5D + 2)y = 260 sin 2t. 



(a) Find a particular solution. (Answer. 11 cos 2t — 3 sin 2t.) 

(b) Show that the corresponding characteristic polynomial is factorized as 

p{X) = (A2 + 3A + 2)(a2 + A + 1), 
and hence the zeros have negative real parts. 

(c) Show that the determinant 



satisfies the Routh-Hurwitz criterion. 
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